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I. INTRODUCTION 



M. K. VolkovE E. A. Kuraevjt] and Yu. M. Bystritskij^ 

Joint Institute for Nuclear Research, Dubna, Russia 

The two-photon decay widths of scalar mesons <r(600), /o(980) and do (980) as well ao — > p(o>)7 
and /o — > /o(w)7 are calculated in the framework of the local Nambu-Jona-Lasinio model. The 
contributions of the quark loops (Hartree-Fock approximation) and the meson loops (next 1/N C - 
approximation where N c is the number of colors) are taken into account. These contributions, as we 
show, are the values of the same order of magnitude. For the fo decay the K-\oop contribution turns 
out to play the dominant role. The results for two-gamma decays are in satisfactory agreement with 
modern experimental data. The predictions for ao — > p{^)"i and fo — > p(w)7 widths are given. 
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\Q ' In recent paper the radiative decay <j> — > /07 width within the local Nambu-Jona-Lasinio (NJL) model 0, 0, 
EL S B S B B E3] has been calculated. In this work, same as in we took into account not only the quark loop 
contributions (Hartree-Fock approximation) but also the meson loop contributions (next 1 /-/V c -approximation where 
N c is the number of colors). We should note that meson loops give the contribution of the same order as the Hartree- 
O ■ Fock approximation due to fractional charge of quarks and integer charge of mesons. The gauge invariance of the 
amplitudes leads to absence of ultraviolet divergences in the relevant loop integrals. Thus, the explicit dependence of 
these amplitudes of the external momenta was obtained [Tlj |. 

In this paper, we will consider the two-photon decays of the scalar mesons cr(600), /q(980) and a (980) as well as 
radiative decays ao — + p(y)l and fo — > p(u})j- 

In the case of the quark loop we consider only a real part of the relevant loop integral. This prescription permits 
us to take into account the condition of the "naive" quark confinement. Some theoretical arguments supporting this 
procedure can be found in [l2T |. As for the meson loops, both the real and the imaginary parts were taken into account. 
OO I The structure of our paper is the following. In Section [Til the NJL quark-meson Lagrangians, corresponding 
. parameters and the coupling constants of our model are defined. In Section [TTTT the methods of quark and meson loop 
calculation are described; the contributions of quark and meson loops to the amplitudes and the widths of two-photon 
^-jl . decays of the scalar meson are presented. 

' In Section HVl we consider the ao — > p(oj)'y and fa — > p{oS)"f decay widths. 
0^ ! In Section [Vj we discuss the results obtained. 

o 



II. LAGRANGIAN OF THE NJL MODEL 

The initial four-quark Lagrangian of the local Nambu-Jona-Lasinio model has the form 0, @, : 



£ = q id + eQA-Mo 



q 



+ % + (9(*A a 75)g) 2 ) - ^ ((f7„ A Q g) 2 + (q( W \ a )q) 2 ) , (1) 

where q = (u, d, s) , and u, d, s are the quark fields, Mq = diag (to°, m° d , is the matrix of quarks current masses, 
Q = diag (2/3, —1/3, —1/3) is the quark electric charge matrix, e is the elementary electric charge (e 2 /47r = a = 
1/137), Aj, i = 1 . . .8 are the well-known Gell-Mann matrices and Ao = \/2/3 diag (1, 1, 1), G and Gy arc effective 
coupling constants of four-fermion interactions. 
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The procedure of bosonization and the renormalization of meson fields leads to the effective quark- meson Lagrangian. 
The part of this Lagrangian which we will use is the following 0, 0, HI : 



C 



eff 



id + eQA - M 



q + q 



g<j u K<?u + ga a X s a s + ga u X 3 a + 1755^ (K+n + + X^-tt ) 
9p 



+ H59K {^k+K + + X K -K ) + y (X 3 p + X u uj) 



q, 



(2) 



where M = diag (m u , m^, m s ) is the constituent quark mass matrix and m u = = 263 MeV, m s — 406 MeV, 
X u = (V2X Q + A 8 ) /V3, A s = (-A + \/2A 8 ) /-y/3, X„± = (X 1 ± iA 2 )/V2 Xk± = (Xa ± «A 5 )/V2. Taking into account 
the six-quark interaction of t'Hooft leads to mixing of a u and a s states [7l. la. [To| . Scalar isoscalar mesons /q, <j are 
the mixed states 



fa = <j u sin a + <j s cos a, 
a = a u cos a — a s sin a, 



(3) 



with the mixing angle a = 11.3° 

The coupling constants from the Lagrangian ([2]) are defined in the following way 0: 



9^ 
9<y s 

.971 

9K 
9,> 



(4/ A (m u ,m u )) 
(4/ A (m.,m.))' 
= 2.84, 

F 

1 TV 



-1/2 



-1/2 



= 2.43, 
2.99, 



m u + m s 



= 3.01, 



2F K 
V6g au = 5.95, 



where we use the Goldbcrgcr-Treiman relation for g n and gx constants, F v = 92.5 MeV and Fk = 1.2 F^, and 
/ A (m, m) is the logarithmically divergent integral which has the form: 



I(m, m) 



(fc 2 + m 2 ) 2 



In 



A 2 



A 2 



N c = 3. 



(An) z V V m ^ 7 A 2 +m 2 /' 
This integral is written in the Euclidean space. The cut-off parameter A = 1.27 GeV is taken from [3, 



III. TWO-GAMMA DECAYS OF a(600), /o(980) AND a„(980) 



The amplitudes of the 27 decay can be expressed in terms of the quark and meson loop integrals. 

The quark loop contribution to the amplitude is given by two triangle type Feynman diagrams (see Fig. [TJ a): 



rpq 



+ ((tfijM) ^ (92,^)) • 
Applying the Feynman procedure of joining of the denominators 

1 



d 4 k s p 


7„ (k + m q 


7m 


[k + qi + 


m q) 


(k - 92 + m q 


)] 


™ 2 (fc 2 - m 2 ) ((fc + ft) 2 - to 2 ) 


((* 


- 92) 2 - ™ 2 ) 





(4) 



(fc 2 - m 2 ) Uk + qif - mfj ((k - q 2 f - mfj 



1 



1 1 

= I dx I 2ydy 

{ { {{k-q x y?-(ml + qW)) 



(5) 
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a b c 

FIG. 1: The Feynman diagrams of quark and meson contributions to two-photon decays of scalar mesons: S —* 27. 



where q x = xq 2 — xq\, x — 1 — x, we obtain for : 



T L = - (g»v (3192) - qiv<i2n) T q , 

TT 



T q = 2m q dx ydy 



1 — Ay 2 xx 
m 2 — M^y 2 xx 



(6) 
(7) 



For meson loops an additional Feynman diagram with two photon- two meson vertex contributes as well (see Fig. [2 b, 
c). To restore the general gauge invariant form of the amplitude, we can nevertheless consider only two triangle type 
Feynman diagrams: 



(2k + gi ) (2k - q 2 ) v 



a/;M - . .. 

(k 2 - M 2 ) Uk + qif - M 2 \ l(h - q 2 f - M 2 
+ Hqi, M) (Q2,v)) ■ 
Extracting the term ~ q\ v qi^ and adding the relevant term ~ we obtain: 



(8) 



with 



T M = 2 I dx I ydy 



y 2 xx 



Standard evaluation of these integrals leads to 







jiq 



M 2 - M 2 s y 2 



-H4 



Sh 



T 



.,M 



4M 2 1 s J ' 



where z% = 4m 2 q /M 2 , zf = 4M 2 /M 2 , 



F(z) = Re[l + (1 - z)$(z) 
$(*) = z<f>(z)-l, 



r 2 _ ln 2 



2 



4- In 1+v lZf 



( arctan ■ 



z < 1, 



z > 1. 



(9) 
(10) 

(11) 
(12) 



(13) 



We remind that for the quark loop contribution the imaginary part of the function $(z) must be omitted (see [i2j]). 
and for the meson loop contribution both the real and possible imaginary parts are relevant. 
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Similar expressions were obtained in [11] , where the imaginary part of the quark loop contribution was taken into 
account. 

The vertices of the quark- meson and quark-photon interactions were given in ([2]). The vertices of the meson-meson 
interaction in the framework of the NJL model have the form (for details see [HH, 0|) : 

V^k+k- = V aoK+K - =-2(2m u -m s )^, (14) 

9<Tu 

V asK +K- = 2V2(2m s ~m u )^, 

9<7s 

g 2 

V a ^+w- = -4m„ — . 

9<Tu 

The general structure of the two-photon scalar meson decay amplitudes has the form 

7s 77 = — {g^u (9192) - qiq%) as 7 7- ( 15 ) 

7TTYlq< 



The expression for the width has the form: 



Ml a 2 g 2 a 2 

^ = 64^« losrrl ■ (16) 



A. The decay ao — > 2j 

The amplitude a aa ~ n of ao — > 77 contains the contribution of u, d quarks and the -ftT-meson intermediate states. 
The color-charge factor associated with u, d quarks is N c (| — |) = 1. Thus, 

<f 77 = (17) 
Taking the Zf-meson loop contribution we obtain: 

«a o77 — a a o77 "+" a a o77 — 

{ ao) g« u ±g* u M 2 K ao { ao) 

= 0.482 - 0.114 = 0.367. (18) 



The corresponding width is 

r ao (980)^ 77 = (2.25 keV)|a Q077 | 2 = 0.29 keV. 
This value is in satisfactory agreement with the experimental values (see Table UJ. 



B. The decay / (980) -> 77 

In the case of the / — > 77 decay we also have the contribution of u, d and s quarks and the _RT-meson and pion 
intermediate states. We should recall that /o(980) meson consists of two components: a u and a s (see ©). The color- 
charge factor associated with u, d quarks is N c (f + 5) = § f° r Cu-component of /o and N c ^ = ^ for cr s -component 
of /o- Taking the K -meson and the tt- meson loop contribution we obtain 

a/0-77 = 3^ [ZfJ Bin a - —F (z fo ) [^—^ J cos a + 

+ g^t^ 2V ~ H2ms ~ m " )c ° sa ) ^ {zl) ~ 
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2 2 

m U 9tT 7T ^ ( IT \ 

sma——z fo <P(z fo ) = 



For the width we have 



= 0.157 - 0.417 - 0.022 + 0.589 + 0.082 - 0.038i = 

= 0.385- 0.038i. (19) 



r M 98o)^rr = (2.25 keV) |a /o77 | 2 = 0.33 keV. 



This value is also in satisfactory agreement with the experimental values (see Table [ 

C. The decay a — » 77 

In the case of the a — > 77 decay we have 

' 9a B m u 



- \ F ( z a) cos a + (z s a ) [ t " 7s "" u ) sin a 

( 9k m u „ /„ x 

? 777T 2 ( 2m « ~ TO s) cosa + 



+ 9K JT72 2 ^ (2m s - m u ) sin a ) z* $ (zf ) - 



-,2 

= 1.89 + 0.057 - 0.041 - 0.043 + 0.92 - 0.98i = 

= 2.78- 0.98i. (20) 

Let us notice that unlike /o(980) — > 27 decay, where contributions of u and d quarks have opposite signs, here they 
have the same sign. 

The corresponding width is 

r ff(6 oo)-. 77 = (0.51 keV)K 77 | 2 = 4.3 keV. (21) 

The experimental value of the mass and the width of the a meson is not well established. We present the width of 
cr for two other masses: M a = 450 MeV and M a = 550 MeV. They are 

Ta(450)^77 = 2.18 keV, 
r CT (550)-*77 = 3.53 keV. 

The comparison of our results with the experimental data and some other model predictions is given in Table [T] 

IV. THE DECAYS a -> u>(p)j, fo -> w(p)7 

In this section, we will consider the following decays [4(| : 

a Q (p) -> oj(q) +7(fci), 
ao(p) -> p(q) +l{ki), 
fo(p) -> v{q) +j(ki), 

Hp) -> P (q) + i(k 1 ), 

P 2 = M 2 S , q 2 = M*, fc 2 = 0, (22) 



where Ms = Mf 0tao = 980 MeV is the mass of decaying scalar meson [13| and My = M u , p is the mass of vector 
meson. Now the matrix element in the general case has the form: 

Mi = e (<7^ fc iM - g^v (qki)) e 7 e y , (23) 

e 7 = e(fci), e v = e(q), (24) 
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V(q 2 ) 



FIG. 2: The Feynman diagrams of quark and meson contributions to decays of scalar mesons to photon and vector meson: 

where i = {ao — > W7, ao — > fry, fo — > W7, fo —> pj}, and Ai contains all vertex constants and dynamic information of 
the processes. The radiative decay width then has the form: 



32MJ 



(25) 



To calculate quark and meson loop contributions to coefficients A4, we need quark- meson (sec J2])) and meson- meson 
vertices (see (fT4|) ). and vertices of interaction of vector mesons with the pseudoscalar ones (3. M. Il4|: 



9ui>K+K- — 9p»K+K- 



9 P 



0+ 



g P ^7T+7T- = 9 P ip+ -p-Y ■ 

A. The decays ao —* uj"f, ao —* p7 



Let us consider the decay of the isoscalar meson ao with w-meson production. The amplitude of this decay will 
contain contributions from the quark loop and if-meson loop. The 7r-meson loop is forbidden since the ao — > 2-7T 
vertex is absent. 

Quark contribution consists of u and d quark loops (see Fig. [2l a) (since u> does not contain s quarks) : 



9p_ r (uA) 





^ + k + m^j 


\k-k\ + m u 


e 7 (k + m^j 




27T 2 | 


(q + kf - ml 


) (k 2 - ml) ((fc - fci) 2 - ml 


) 



where C^'Su^ — 3g CTu , where factor 3 = iV c is the color factor. 

Standard Feynman procedure of denominators joining and loop momenta integration leads to: 



M^U = e | Ctdl 7 Re (I u ) (q v k llt - g„ v (gh)) %e£, 



where 



1 1 

4m u dx dy- 
J J 1 





y (l - Ay 2 x (1 - x)) 



i\ — y(l — y)(l — x)q 2 — x(l — x)y 2 p 2 + it 



(26) 



(27) 



(28) 



Let us consider now kaon loop contributions. The kaon contribution consists of diagrams of two types (see Fig. [2]b, 
c). This contribution can be written in a form similar to (1241): 



(29) 



where C, 



,(K) _ 

a ^Ljj — 9aoK+K- 



and 



I K I dx I dy ^y 2 x{\ - x)y 

Mjj- — y(l — y)(l — x)q 2 — x(l — x)y 2 p 2 + ie 



(30) 
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More details of this type of diagram calculation can be found in [l[ . 
The amplitude of the process ao — ► tv"f then has the form 

Map-H^ = e A ao ^^ {qvkip - g^ (qh)) e^e%, (31) 

where 

A ao ^ = 3g au Re(I u ) + g aoK+K -I K = -1.78374+ 0.159415 = -1.62433. (32) 



The decay width is: 



T ao ^ 7 = 115 keV. 



The decay ao — > pj can be considered in complete analogy with the decay clq — > ary. Let us note that the quark 
contribution to a — ► fry is three times smaller than the quark contribution to the ao — > cvy decay. However, the kaon 
loop contributions are the same. As a result, the amplitude of the ao — > p"f decay has the form: 



2 

A ao ^ P1 = .ga u Re (/„) + g aoK +K-lK = -0.598209 + 0.156921. 



Afao-^ = e -^A ao ^ pl (q„fci M - g^ {qh)) e"e%, (33) 



The decay width is 



r Q0 ^ P7 = 8.5 keV. 



Let us note that in both these processes the main contribution comes from quark loops. 

B. The decays fo — * luj and fo — > pj 

The total amplitude for the fo — > uj decay has the form: 

A h ^ UJ = sina A^ U7 +sina Af u ^ W7 + cosa Af a ^ W7 = -0.115 + 0.0195+ (-1.068) = -1.1635. (34) 
The decay width is then: 

r /o ^ 7 = 60keV. 

Let us consider the decay fo — > fry. The quark loop contribution to a u component here exceeds the relevant 
component of the /o — > decay: 

A f -*m = sin « ^"'5 P7 + sin « A l u ^fr, + sin « A a u ^ P1 + cos a A*^ = 

= -0.353 +(0.0436 -0.00542*') + 0.027 +(-1.05) = -1.3324 -0.00542i (35) 

The fo — ► fry decay width is then: 

r /o _ P7 = 82keV. 

Note that in both the cases (the decays of /o(980) to /try or ujj) the main contribution comes from the kaon loops 
related with the a s component. 

Unfortunately, at present we do not have any experimental data for this decays. 

V. CONCLUSION 

The calculations of the radiative decays in the NJL model show an important role of both the quark and meson 
loops. Moreover, for the fo meson decay the kaon loop turns out to provide the dominant contribution. This can be 
understood if one takes into account that in the 27 decays the fractional quark charge gives effectively factor 1/N% 
into quark loop contribution. It is worth noticing that the situation here is similar to the one that takes place in the 
case of (ft — > /07 decay [l[ and ao — > p(w)y and fo — > p(oj)j. This fact permits one to understand the success of such 



17 and , 

models as the model of a kaon molecule [IH, [lj| as well as the four-quark model [171 11 
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The description of radiative scalar isoscalar meson /o decays was also considered in the linear u-model [HI, H(| [2l[ • 
However, in these papers the scalar meson ao(980) and a radiative decays were not considered. 

The NJL model used here allows us to take into account both the quark-antiquark state which manifests itself in 
the form of quark loops, and the hidden four-quark state which shows up as meson loops. Let us emphasize that in 
the framework of the standard NJL model we can describe the radiative decays without any additional parameters. 

The results obtain in this paper (see Table HI arc in agreement with the existing experimental data. Unfortunately, 
the data for the do — > p{oS)^f and fo — > p(w)j decays are absent up to now and our results for these decays can be 
considered as a prediction. 
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TABLE I: The table of two-gamma decays of the scalar mesons <r(600), /o(980) and ao 



r (a -* 77) {exp.), keV 


r (ao — » 77) (theor.), keV 


0.30 ± 0.10 [22] 


0.29 (This paper) 


1 (jo — > 77; (exp.), Kev 


1 ^jo > 77J (ineor.), Kev 


0.42 (Solution A) [23] 


0.33 (This paper) 


0.10 (Solution B) [23J 


0.21 — 0.26 [16J 


0-2051^^1^;^ [24] 


0.22 [25J 


n oq+0.09 rofil 

0.28_ 13 126] 


0.33 I2YI 


42 ± 06 ± 18 [281 


31 [291 


0.29 ± 0.07 ± 0.12 [30] 


0.28+S ?? [311 

-O.Li 1 


0.31 ± 0.14 ± 0.09 [32] 


0.20 [33] 


0.63 ± 0.14 [34] 


0.24 [35] 




0.27 [36] 


1 (<r — ► 77j (c^p.), keV 




3.1 ± 0.5 (Solution A) [23| 


4.30 (M CT = 600 MeV)(This paper) 


2.4 ± 0.4 (Solution B) [23] 


3.53 {M a = 550 MeV)(This paper) 


4.1 ±0.3 [37] 


2.18 (M a = 450 MeV)(This paper) 


3.8 ±1.5 [26] 




5.4 ±2.3 [34] 




10 ±6 [38] 





